Magnetoresistance of a Wigner liquid in a parallel magnetic field 
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It is assumed that in a two-dimensional electron system with strong correlation (a Wigner liquid), 
appearance of some relatively slow-moving objects (clusters) composed of small number of electrons 
is possible. Such clusters may exist in addition to ordinary mobile carriers of the Fermi type. 
They can pin to inhomogeneities and play the role of additional scatterers. The clusters composed 
of two and three electrons are discussed (for near order as in triangular lattice). The number of 
the clusters depends on temperature and a parallel magnetic field (so accordingly for conductivity 
and magnetization). In the frame of a simple model, a resistivity increasing and a metal-insulator 
transition with an increasing of a magnetic field are proved. Near this transition, the resistivity 
changes with temperature according to a linear law. The model gives a nonlinear dependence of 
magnetization on a magnetic field. 



PACS: 71.27.+a 

In recent years, low-density systems of two-dimensional 
electrons (holes) are intensively studied (see, e.g., 
reviews 1,2 ). That sort of systems exhibits a number 
of unusual properties. One of the main properties is a 
metal-insulator transition which occurs upon a decreas- 
ing of the carrier density. Moreover, with an increasing 
of a parallel magnetic field, a growth of resistivity and a 
metal-insulator transition are observed. The present pa- 
per is devoted to the theory of such systems in a parallel 
magnetic field. 

The above-mentioned phenomena may have various 
reasons. For example, theoretically the effective mass 
m* could increase with spin polarization, but the experi- 
ments apparently do not manifest that. Further, interac- 
tion of carriers with impurities can increase with spin po- 
larization. This effect was considered in the work 3 - where 
the interaction with charged impurities was discussed. 
However, the concrete realization of the idea in this work 
gives rise to doubts. The point is that, for the screened 
interaction with charged impurity, the used in the work 3 
formula is applicable only in the limit r$ << 1, while, on 
the contrary, the limit r$ >> 1 takes place for the con- 
sidered systems (rs is the dimensionless average distance 
between the particles). In the worki, another mecha- 
nism connected with the presence of hole traps near the 
SiC>2/Si interface and their possible recharging was dis- 
cussed. Though this mechanism is possible but its effect 
depends on a sample and changes from specimen to spec- 
imen, whereas experiments show more or less universal 
behavior for any utilized sample. Seemingly, it is neces- 
sary to find some universal mechanism. 

Apparently, the universal reason exists because of some 
universal property of considered systems. Namely, a 
low-density system of two-dimensional electrons (under 
rs >> 1) is the system with the strong correlations, 
i.e. it has a near order in particle arrangements (as in 
a Wigner crystal and, therefore, such a system is called 
a Wigner liquid sometimes). The use of this property is 
the main point of the present work. 

First of all, it is necessary to do some remarks. In the 
limit rs >> 1, there is the following hierarchy of ener- 
gies. 1) The Coulomb energy is the greatest one and is 



taken as the unity. After all, it is this energy that results 
in a Wigner crystal and gives a near order in a Wigner 
liquid. 2) The energy of zero oscillations (or equivalently 
of typical plasma oscillations) is of the order of l/y/r~s~. 
Apparently, because of the zero oscillations the liquid do 
not crystallize right up to rs ~ 37 (see the work£). 3) 
At last, the Fermi energy is of the order of 1/rg if this 
energy is estimated by using of band effective mass mj, 
as in a Fermi gas. In reality, the Fermi energy is even 
smaller because the effective mass m* renormalized due 
to interaction is greater than nib- Further, we shall deal 
with the energies of the order of the Fermi energy that 
is smallest one among the all energies and is of the order 
of exchange energy (see below) . 

The main idea of the present work is connected with 
the conception accepted in the work 6 . In this work, an 
estimation of the exchange interaction of particles (and, 
accordingly, of the effective mass m*) was done. For 
that, two neighboring particles interacting with external 
potential produced by environment and with each other 
via Coulomb were considered. Either of the two takes up 
its position mainly in its potential minimum, and from 
time to time they interchange their positions. With the 
help of this model, the exchange splitting of energy levels 
(Ea — Es) (its dependence on rs) has been found. The 
exchange Hamiltonian that gives the splitting may be 
presented in the form: 

H ex = A S |(S 1 S 2 ) + (Q 1 Q 2 ) + 

4(S 1 S 2 )(Q 1 Q 2 ) - -} , (1) 

where the indices 1,2 are particle numbers, S is the spin 
operator, and Q is the quasi-spin operator that corre- 
sponds to the two- valley case (Q = 1/2). Further, the 
two- valley case as at Si02/Si inversion layer is exam- 
ined. The constant (-1/4) in the brackets is added for 
convenience. The value As and the effective mass m* are 
connected with the energy splitting in the following way: 

n 

2As = Ea — Es ~ — - , (2) 
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where n is the carrier concentration. The equation (J2J 
gives estimation of m* . Note that here the right-hand 
part is of the order of the Fermi energy of quasi-particles. 

Starting idea. Further, one would think, it is nec- 
essary to operate in accordance with the Landau theory 
of the Fermi liquid (by using the effective mass estima- 
tion J2J). However, the question arises about adequacy of 
similar approach in our case. Actually, one may imagine 
a situation when, for example, a pair energy is smaller 
than the double Fermi energy of quasi-particles, and so a 
number of the pairs would arise together with the quasi- 
particles. 

In order to elucidate the situation, let us turn to the 
physical picture which was implied in the work&. In this 
work, the exchange energy of two nearest neighbors was 
calculated. If take into account all jumps of a particle in 
liquid, then one can see that the minimal energy of a par- 
ticle can be even lower than in a pair, and this decrease is 
of the order of the exchange energy. Obviously, this min- 
imal energy can be attributed to quasi-particle with zero 
momentum. So, as long as the all relevant to our problem 
energies are of the same order, then, in particular, it is 
permissible to compare a pair energy and the Fermi en- 
ergy of quasi-particles. Therefore, the above-mentioned 
effect can exist and then influence of the pairs on the 
properties of the system should be taken into account. 

One can consider instead of the pairs some other for- 
mations. For near order as in a Wigner crystal (i.e. as 
in triangular lattice), a cluster composed of three near- 
est neighbors is the most appropriate object. Obviously, 
energy of such a cluster (per a particle) is smaller than 
energy of a pair, and so, probably, it is necessary to do a 
preference for a three-particle cluster. 

It is unlikely that one should consider more compli- 
cated objects in liquid because for that the ordered ar- 
rangement of particles is necessary not only for nearest 
neighbors but for next nearest neighbors as well. So fur- 
ther we shall consider the pairs and the threes only. 

Evidently, the number of the clusters depends on a 
magnetic field (since clusters have spin) and tempera- 
ture and the same relates to the number of mobile car- 
riers of the Fermi type (fermions). It is natural to con- 
sider that the clusters are slow-moving formations be- 
cause their jumping amplitude is smaller than for a par- 
ticle. It means that they can be pinned even by small in- 
homogeneities and be excluded from the electric current. 
Moreover, there is an additional contribution to a resis- 
tance due to fermion scattering by the pinned clusters. 
Therefore, the resistivity of the system should depend on 
a magnetic field and temperature as well. The above pic- 
ture is a basis for proposed description of the properties 
of the system. 

In some extent, our approach is partly similar to the 
model used in the work 4 . However, unlike this work, we 
appeal directly to the intrinsic properties of the strong 
correlated system rather than to the external factors. 

Clusters. First of all, it is necessary to describe the 
properties of the clusters. The states and energies of the 



clusters are meant. Let 5* and Q be the total spin and 
quasi-spin of a cluster, accordingly. We start with a pair. 
The quantum numbers corresponding to its minimal en- 
ergy are S = 1, Q = or S = 0, Q = 1. With the help 
of ||TJ, it is easy to calculate the minimal energy E% of a 
pair with a result: 



E 2 /2 = -3A s /8 



(3) 



(here the energy per a particle is indicated). 

Now consider a cluster of three particles. The state 
with S — 3/2, Q = 3/2 has the maximum energy and 
is of no interest for us. Passing to consideration of the 
states with 5" = 3/2, Q = 1/2 or S = 1/2, Q = 3/2 with 
the same energy E' 3 , we obtain: 



E' 3 /3 = -7A S /12 



(4) 



(per a particle as well). These states are doubly degen- 
erated for the fixed quantum numbers S, S3, Q, Q3. At 
last, for the energy E' 3 ' of a state with S = 1/2, Q =1/2 
we have: 



E' 3 '/3 = -3A s /4 



(5) 



This state is not degenerated for the fixed quantum num- 
bers S, S3, Q, Q 3 . 

Model. Here we restrict ourselves to the simplest 
model. For the determination of the equilibrium prop- 
erties, pinned clusters and a gas of mobile carriers 
(fermions) are considered without any interactions. An 
interaction of fermions with clusters is taken into account 
only for estimation of resistivity. The energy of the sys- 
tem is written as: 



Ha 



»v(p) 



p,cr 



£ [7(£)£ Q (£, v) + HS 3 ] iV(£, v) . (6) 
£,1/ 

Here the first sum relates to the fermions (the summa- 
tions take places over the momentum p and spin projec- 
tion a = ±1/2). The second sum relates to the clusters. 
The following notation for the set of quantum numbers 
is used: 

(5,5 3 ,Q,Q 3 )^S . 

In the second sum, the summations occur over these 
quantum numbers and state numbers v. The value 
7(E)i?o(S, v) is a cluster energy (with the quantum num- 
bers E and in the state v) and 7(E) is number of particles 
in the cluster. The constant takes into account ad- 
ditional degeneracy of a level (k = 2 for the states with 
the energy (0J and n = 1 for all other states). Finally, 

H = g*fiBB . 

Here B is the magnetic field, g* is the effective g - factor 
(supposed to be the same for the fermions and clusters), 
fiB is the Bohr magneton. 
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Further it is meant that only one cluster can be lo- 
cated in the every state v. Apparently, the total number 
of these states is of the order of the total number of the 
particles. One would think, that in a liquid it is neces- 
sary to consider a cluster as a particle with enough big 
effective mass and characterize it by momentum. How- 
ever, in our opinion, it is not unreasonable to describe 
the clusters by localized states because of imperfections. 

Equilibrium properties. The cluster characteristics 
can be obtained with the help of statistical sum Z(y) for 
a given state: 



exp 



Z{y) - 1 

7 (E)[£b(S,i/)-A*] +HSs 



T 



(7) 



where /i is the chemical potential of the system. 

The fcrmion distribution function has the usual form: 



r(p) = { 



e(p) - 
T 



l} 1 , n„ = n-Ha. (8) 



For quadratic dependence of the fermion energy on the 
momentum (e(p) = p 2 /(2m*)), we have the following 
expression for the fermion concentration: 

n a = I^n CT (p) = ^ln{l + exp(/i CT /T)} (9) 



(for a fixed spin projection and for a fixed valley). 

As usual, the chemical potential \x is defined by the 
condition that the total number of the particles is equal 
to the sum of the numbers of fermions and electrons in- 
side the clusters. In the case of two valleys, this condition 
looks as: 

n = np + nr , np = 2 n„ , 

cr 

-T 5Z(v) 



nr — 



E 



Z{y) 8E Q {V,v) 



(10) 



Here n is the total concentration of electrons, rip is the 
concentration of mobile carriers (fermions), and nr is the 
concentration of electrons inside the clusters (it is deter- 
mined by a functional derivative of the function Z(y), see 
the expression J7J|). 

Up to now, the various expressions were written in gen- 
eral form that took into account the clusters containing 
the different numbers of electrons. Further, we restrict 
ourselves to clusters with three electrons for which the 
more low energies are obtained than for pairs (sec JSJ - 
|J5J). For these clusters, it is necessary to take into ac- 
count not only the states with lowest energy (JjJJ but also 
the states with energies (J3J because, in a sufficiently high 
magnetic field, a contribution of clusters with spin 3/2 
is able the main. Furthermore, for simplicity, we neglect 
with dispersion of the cluster levels, i.e. we consider that 
every energy Eq does not depend on the state number v. 



Instead of the designation -E'o(E) for which it is necessary 
all the time to indicate a set of the quantum numbers, 
further it is convenient to use some other designations, 
namely: 



E (1/2, 1/2) -»e , 
£ (l/2,3/2)=£ (3/2,l/2) 



fll) 



(the arguments correspond to the spin and quasi-spin 
values). As appears from the expressions 10} and JSJ, 
e < e x . 

After these remarks, the statistical sum (JJJ is written 
in the form: 



8 exp 



Z = 1 + 4 exp 



-3(ei - /i)/T 



3(e - n)/T 
3H 



cosh • 



cosh^ 
H 



2T 



3 cosh ■ 



2T 



(12) 



For the concentration of electrons inside the clusters, we 
have: 



nr = 3no- 



1 



Z 



(13) 



where no is the concentration of the cluster states. 

Now let us find the spin polarization M depending on a 
magnetic field. Its ratio to the maximum possible value, 
M m = n/2, can be written in the form: 




M 

Mm " 

cosh [(/i 



2T 



V 



dZ 
~dH 



3Z 

ff/2)/(2T)] 



cosh[(/x-F/2)/(2T)] 



(14) 



Here the first term in the right hand side is the cluster 
contribution, and the second term is the fermion con- 
tribution. The constant n corresponds to the maximum 
part of the electrons which can be located inside the clus- 
ter states (apparently, one can think that n is less and of 
the order of unity) , ep is the Fermi energy for the two- 
valley case at T — and without taking into account any 
clusters: 



3n n 



V = 



2m* 



(15) 



Zero temperature. Let eo < tp. In a magnetic field 
and at T = 0, it is necessary to take into accounts only 
the lowest levels of the clusters, namely: 



e - (H/2)/3 , ei - ei - H/2 . 



The energy shifts for electrons are 5e^ \ = ±H/2. As a 
result, we have the following expressions for concentra- 
tions of fermions with the different spin projections and 
electrons inside the clusters (with minimal energy): 



n } = —(e -2H/3) , 

7T 
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in 

H = —(eo + H/3) , 

7T 

777 

n T = n (2e -ff/3) ; 

7T 

(0 < H < 3e /2) . 
In the greater magnetic fields, 

3e /2 < H < 3(2e F - e ) , 



(16) 



(17) 



a devastation up to zero of the cluster levels and an in- 
crease of nj — > n take places (according to the linear 
laws). These behaviors are without taking into account 
the levels with the energy t\ that is in the case t\ > 2e F . 

Now, let e < £i < 2e F . The level with energy e x comes 
into the play at magnetic field H/3 > t\ — eo- There are 
the two cases. Firstly, this level can be connected up 
under nj — 0, i.e. when the magnetic field falls into the 
interval i|17|) . In this case, the cluster levels are devas- 
tated not completely (at the left border of the interval 
Q17JI. i.e. under t\ = 3e/2, these levels are not devastated 
at all). Secondly, this level comes into the play under 
n-f > 0, i.e. under ei < 3e /2. Let us discuss the latter 
case. 

In the range of magnetic field < H < 3(ei — eo), there 
are the previous expressions (|16f) . In the point H/3 = 
ei — 6q, the states with initial energy ei are filled up 
instead of the states eo and, without taking into account 
a dispersion of the cluster levels, that takes place by a 
jump. In the range of magnetic field 3(ei — eo) < H < e±, 
we have: 



n T = (ei - H) 

7T 



-ei 



(18) 



n r 



(2ci - H) 



That continues up to the magnetic field H = ei, i.e. up 
to the complete polarization, when the value of n-\ equals 
to zero and the values of n^ and nr become saturated. 
The dependence of any quantity, including a resistivity, 
on a magnetic field has a form of a broken curve, and for 
a magnetic moment even with jumps. It is clear that the 
curves are smoothed by temperature. 

We note that the above expressions are valid until the 
obtained number of electrons inside the clusters is smaller 
than a maximal possible one, i.e., until nr < 3no- That 
takes place under the condition 

, ei 

1 -2e7 < ^ 

Now about the magnetic moment. Let us consider the 
case eo < e F , eo < ei < 3eo/2 (see ijTflfr. (TT%1) for the 
two intervals of a magnetic field). We obtain for the low 
magnetic fields: 



M _ 5H 1 / e 
M m 9e F 3\ e F 

< H < 3(ei - e ) 



(19) 



For the high magnetic fields up to the complete polariza- 
tion, we have: 



M , H-ei 
3(ei - e ) < H < e x . 



(20) 



The jumps of the magnetic moment take places at the 
points H = and H — 3(ei — eo). The ratio of suscepti- 
bilities in the two ranges of low and high magnetic field 
is 



Xi 

X2 



In considered case, the complete polarization takes place 
under H = e\ < 3eo/2 < 3e F /2, while without the clus- 
ters under H — 2e F - 

In connection with the jump of magnetic moment near 
the H = 0, it is necessary to note the following. In 
the workl, the direct measurements of the thermody- 
namic spin susceptibility were produced. It is revealed 
that magnetization behaves by nonlinear manner, the 
spin susceptibility becoming greater under magnetic field 
lowering. Possibly, that is the effect of the clusters. 

Resistivity. For the resistivity p, we use the ordinary 
expression: 

(21) 



n F e 2 r 



where r is the relaxation time. The question is what 
is necessary to take as the r? If an interaction of the 
mobile carriers with impurities is weak (not so with the 
clusters) or the number of the impurities is sufficiently 
small, then the relaxation time is proportional to an im- 
purity concentration, i.e., for the scattering by clusters, 
1/r ~ nr. That is the main for the sufficiently pure spec- 
imens. Perhaps, it is necessary to take into account the 
other scatterers too, then, in the simplest case, one can 
write: 



1/r ~ nr + a n», 



(22) 



where n, is the concentration of the other scatterers and 
the factor a takes into account a difference of the scat- 
terers. As a result, taking into account the expressions 
(|13|l and (|1U|) . one can write: 



r)(Z-l)+AZ 

(i - V )z + v 



A 



(23) 



Let us consider the low temperatures T << e F . Under 
H = 0, from the expressions (|10|) and l|13|l . for determi- 
nation of the chemical potential, we have the equation: 



2m* fi 



3n 



(1/4) exp 3(e - u)/T + 1 



(24) 



(excluding the cluster states with the energy ei > eo). If 
one looks for the chemical potential in the form 

u = e + ui , 
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then one can write the equation (|24[1 as follows: 



V 



£0 = Mi 

e F e F (1/4) exp(-3/xi/T) + 1 



(25) 



Under a partial filling of the cluster levels (for that it 
is necessary eo < £ f) and at T = 0, [i = cq. And that 
is right also at T / for certain values of parameters, 
namely: 



1 = — , Mi = o 

e F 5 



(26) 



If that is the case then the cluster number is a constant. 
In this case, a resistivity does not change while the ex- 
pression (|22|l is valid. So, we have the condition for a 
separatrix. 

For small changes of the parameters, one can ignore 
the first term in the right-hand side of the l|25l) . so that 
one obtains approximately: 



3 '" 1 ~ 1 A 



(27) 



where, for convenience, the designations (3l and (3r are 
used for the left-hand and right-hand sides of the 1)2(1 [I. 
accordingly. From the (|27|l . one can see that, under the 
condition /3r/ '/3l > 1 (mi < 0), the fermion number de- 
creases and the number of scatterers increases, therefore 
the resistivity increases with temperature. Such behavior 
corresponds to a metallic phase. Otherwise, we have an 
insulator phase. In the both cases, the resistivity behaves 
in accordance with a linear law. 

Resistance behavior with temperature depends on a 
value of a magnetic field. Let us consider a magnetic 
field in the interval given for the expression (|18|) . In this 
case, it is enough to take into account the cluster state 
only with the spin 3/2 and quasi-spin 1/2. The previous 
expressions <|24|) - Q26|l. after the replacement 

e -> (e 1 - H/2) , 

remain valid. For example, instead of the l)26|l we come 
to: 



1 



£l - H s /2 Ar, 



Mi = , 



(28) 



where Hs is the magnetic field corresponding to a sep- 
aratrix. In view of Hs < e\, we obtain the following 
condition for origin of separatrix: 

2e F 5 



Under deviation from the l|28(l . the expression for /ii 
coincides with the l)27|) but, in place of the previous left- 
hand and right-hand parts, it is necessary to substitute 
the corresponding parts of the l|28|l . With magnetic field, 
the left-hand part grows whereas the right-hand part is 
constant. Then, as follows from the (|27|l . for the lower 
magnetic field, H < Hs, we have the metallic phase (the 
resistivity grows with temperature) and for the greater 
magnetic field, H > Hs, we have the insulator phase (the 
resistivity decreases with temperature increasing). 

Physical meaning of the appearance of separatrix is 
simple. When a filling of the cluster states is small then 
it increases with temperature. Otherwise, for a large 
filling, it decreases with temperature increasing. The 
cluster filling grows with the magnetic field. At some 
magnetic field, an intermediate case can take place and 
the cluster filling is constant. So, a separatrix appears. 
Thus, when one draws a conclusion from the tempera- 
ture dependance of resistivity then the metal - insulator 
transition can happen. Such a behavior of resistivity cor- 
responds qualitatively to the experiments. 

In the figure, the resistivity dependence on magnetic 
field is shown for different temperatures, namely: T / tp — 
0.05 (dots), T/e F = 0.1 (solid line), and T/e F = 0.15 
(dashed line). The following parameters are taken: 
£ o/ £ F = 0.95, ci/ef = 1-1, 77 = 0, and ^4 = 0. The po is 
resistivity under T = H = 0. The point of intersection 
of the curves corresponds to separatrix. 

The conclusion about enhancement of the resistivity 
under influence of the magnetic field is the direct conse- 
quence of our model accounting for the presence of the 
clusters. The validity of the expressions (12111 and i|22|) is 
assumed. That would be the case if the interaction of 
the mobile carriers with the clusters could be considered 
as weak. For strong interaction, the expression l|22|) can 
be used only in the limit of small number of the clus- 
ters, rtr << n, possibly up to rtr ~ n. However, for 
some values of the nr (of the order of n), that simple 
picture can be incorrect because it is necessary to take 
into account the more complicated phenomena such as, 
e.g., localization of the fermions. 
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FIG. 1: Resistivity dependence on magnetic field for different 
temperatures (see text) 



